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Abstract. We present the first detailed joint modelling of both the timing and spectral properties 
during the outburst decay of transient anomalous X-ray pulsars. We consider the two sources XTE 
J1810-197 and CXOU J164710.2-455216, and describe the source decline in the framework of a 
twisted magnetosphere model, using Monte Carlo simulations of magnetospheric scattering and 
mimicking localized heat deposition at the NS surface following the activity. Our results support a 
picture in which a limited portion of the star surface close to one of the magnetic poles is heated at 
the outburst onset. The subsequent evolution is driven both by the cooling/varying size of the heated 
cap and by a progressive untwisting of the magnetosphere. 
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INTRODUCTION 

Anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs) are two classes 
of X-ray sources whose phenomenology may be explained by the presence of a mag- 
netar\ i.e. an ultra-magnetized neutron star (NS) with a magnetic field as high as 
B ~ lO^'^ - 10^^ G. Both AXPs and SGRs exhibit variability in their persistent X-ray 
emission over a timescale of years, characterized by large flux changes (factor 10-100) 
and drastic alterations in the spectral and timing properties. The recently discovered 
Transient AXPs (TAXPs) represent the far end of this behaviour. Both XTE 1810-197 
and CXOU J164710. 2-455216 underwent a sudden outburst during which their X-ray 
flux increased by more than two orders of magnitude and the emission of SGR-like 
bursts was observed. Their post-outburst decline has continued for a few years after the 
event and has been monitored extensively with XMM-Newton. 

These observations provided a large sample of high-quality data, the analysis of 
which offers an unprecedented opportunity to unveil the mechanisms responsible for 
SGRs/ AXPs activity. Here we present a detailed joint modelling of both the timing 
and spectral properties during the outburst decay of these two TAXPs. We describe 
the source decline in the framework of a twisted magnetosphere model, and make use 
of the spectral models presented and applied by [1] (NTZ in the following), and [2], 
who studied radiative transfer in a globally twisted magnetosphere by means of a 3D 
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Monte Carlo code. Each model is characterized by the magneto spheric twist angle A0, 
the electron (constant) bulk velocity j8, and the seed photon temperature kT. The polar 
field is fixed at 5 = 10^"^ G. Besides, in order to mimick localized heat deposition at 
the NS surface following the activity, we generalize the original application by NTZ by 
allowing for a complex thermal surface map. For details on the method of analysis and 
the observational data, we refer to the original paper by [3] and references therein. 

TAXP ANALYSIS 

XTE 1810-197 was serendipitously discovered in 2003 with RXTE, and a comparison 

with archival RXTE data showed it produced an outburst around Nov 2002, followed by a 
monotonic decline of the X-ray flux [4, 5, 6, 7, 8, 9]. In this investigation we use 8 XMM- 
Newton observations of the source, covering a baseline of 4 years between Sep 2003 and 
Sep 2007. As for CXOU J164710.2-455216, we use 6 XMM-Newton observations made 
from Sep 2006 to Aug 2009, across the Nov 2006 outburst of the source [10, 1 1]. 

Our strategy proceeds as follows. We first jointly fit the pulse profiles in three energy 
bands (0.5-10 keV, 0.5-2 keV, 2-10 keV), allowing for the twist angle, electron bulk ve- 
locity and thermal map to vary in time (but requiring that, for a given epoch, they are the 
same in the different energy bands), and imposing that the geometrical viewing angles 
(see Fig.l) must be the same at all epochs. This step provides us with an estimate of the 
source parameters and viewing angles. We then use our Monte Carlo code to compute 
the phase averaged spectra, at the various epochs, for the same sets of parameters and 
compare with the spectral data. The main reason at the basis of this strategy is that we 
found that spectral fitting alone is unable to constrain the two geometrical angles, while 
lightcurve fitting allows for a better control of complex, time-evolving thermal maps. 

For both sources, we find that the best fit is obtained by using a three temperature 
map, made of a hot spot, a warm corona and a cooler region (see Fig. 1). In both cases 
the method was successful and we reached a satisfactory fit of both the pulse profiles and 
spectra (see Figs. 2,3 and the original paper for details and for the values of the fitting 
parameters). 

The evolution of the single spectral components during the outburst decay was nev- 
ertheless quite different for the two TAXPs. For XTE J1810-197 (see Fig. 4, left), we 
found that the hot cap decreased in size and temperature, until it became indistinguish- 
able from the warm corona around March 2006. The warm corona also shrinked, but it 
initially cooled down to ~ 0.3 keV, and then its temperature remained constant. It was 
still visible in our last observation (Sep 07), with a size down to 0.5% of the NS surface. 
The rest of the NS remained at a temperature comparable with that measured in quies- 
cence by ROSAT Awnng the entire evolution, suggesting that the outburst likely involved 
only a fraction of the star surface [see also 9]. The best-fitting viewing geometry is given 
by;^ = 148°,<§=23°. 

The evolution of CXOU J164710.2-455216 was quite different (Fig. 4, right): in this 
case the hot cap decreased in T and size, and the warm corona remained constant at 
T ~ 0.45 keV while it increased in size. Interestingly, the area of the hot-i-warm region is 
constant (~ 30% of the NS surface), while the remaining ~ 70% is cooler, at ~ 0. 15 keV. 
This may suggest that the thermal map, in quiescence, is made of two components: a 




FIGURE 1. A schematic view of the NS. ^, j are the viewing angles, Q. and jj, are the star spin and 
magnetic axis, the dashed line is the line-of-sight. The star surface is divided into three regions: a hot polar 
cap (red), a warm corona (blue) and a colder zone (grey). Reproduced by permission of the AAS from [3]. 





FIGURE 2. Left: Synthetic and observed pulse profiles for the eight XMM-Newton observations of 
XTE J18 10-197 in the total energy band. Solid lines represent the best-fitting model, dots the observed 
lightcurves. Right: Spectral evolution in the same observations. Solid lines represent the model, while 
dotted lines refer to its individual components. Reproduced by permission of the AAS from [3]. 



warm polar region superimposed to the cooler surface. The outburst might have heated 
a portion of the warm cap, producing the hot zone which then cooled off; we note that 
the hot spot disappeared (~ Aug 08) when the pulse profile switched from three-peaked 
to single-peaked. The best-fitting viewing geometry is = 20°, ^ = 80°. 



Conclusions 

To our knowledge this is the first time that a self-consistent spectral and timing anal- 
ysis, based on a realistic modelling of resonant scattering, is carried out for magnetars, 
considering simultaneously a large number of datasets over a baseline of years. The re- 
sults support a picture in which only a limited portion of the magnetosphere was affected 
by the outburst (twist), and provide identification for the location, meaning and evolution 
of the various thermal components. Future developments will require detailed spectral 
calculations in a magnetosphere with a localized twist which decays in time [12] and 
will be matter of future investigations. 
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FIGURE 3. Same as in Fig. 2 for CXOU J164710.2-455216. Reproduced by permission of the AAS 
from [3]. 
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FIGURE 4. Pai-ameters evolution for XTE J1810-197 (left) and CXOU J164710.2-455216 (right). 
From top to bottom: twist angle (A0), bulk velocity (/3) and area of the warm and hot emitting regions. 
Reproduced by permission of the AAS from [3]. 
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